Beyond quantum Fisher information: optimal phase estimation with arbitrary a 

priori knowledge 
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The optimal phase estimation strategy is derived when partial a priori knowledge on the es- 
timated phase is available. The solution is found with the help of the most famous result from 
the entanglement theory: the positive partial transpose criterion. The structure of the optimal 
measurements, estimators and the optimal probe states is analyzed. The paper provides a unified 
framework bridging the gap in the literature on the subject which until now dealt almost exclusively 
with two extreme cases: almost perfect knowledge (local approach based on Fisher information) and 
no a priori knowledge (global approach based on covariant measurements) . Special attention is paid 
to a natural a priori probability distribution arising from a diffusion process. 
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Quantum states may be employed to improve the 
precision of measurements without committing ad- 
ditional resources such as energy or time [1]. A 
paradigmatic model is the estimation of relative 
phase delay between two arms of a Mach-Zehnder 
interferometer [2]. When the interferometer is fed 
with a light pulse prepared in a coherent state, the 
precision of phase estimation scales as Sip oc 1/ yN, 
where N is the mean number of photons in a pulse. 
On the other hand, by preparation of N photons, en- 
tering the two input ports of the interferometer, in 
a appropriate quantum state, a quadratic precision 
enhancement may be achieved leading to Sip oc 1 /N 
referred to as the Heisenberg limit [3] . 

These idealized results need to be contrasted with 
a more realistic ones when environmental noise and 
experimental imperfections are taken into account 
[4, 5]. For optical implementations the most relev- 
ant disruptive factor is the photon loss. For large 
N, the precision of optimal phase estimation ap- 
proximates to [6, 7] Sip 



i/l — rj/^/rjW, where rj is 
the overall power transmission of an interferometer 
(including scattering, reflections, detector efficiencies 
etc.). The Heisenberg scaling is asymptotically lost 
and the quantum precision enhancement amounts 
to the \/l — rj factor. Although asymptotically the 
prospect of quantum enhanced metrology may look 
bleak, still for moderate N the advantage of using 
entangled states may be very significant [5]. 

Interestingly, the above asymptotic formula was 
obtained independently using two different ap- 
proaches. The first, local approach, aims at maxim- 
izing the estimation sensitivity to small phase vari- 
ations around an a priori known value ip = ipo. The 
main theoretical tool is the quantum Fisher inform- 
ation (QFI), Fq, which defines an asymptotically 



achievable quantum Cramer-Rao (CR) bound [8] on 
estimation precision Sip > 1 / *J Fq . The second, 
global approach, assumes no a priori knowledge on 
the phase — a priori probability distribution is uni- 
form over the [—71", tt) region — and makes use of the 
phase shift symmetry to restrict the class of meas- 
urements to the covariant ones [9]. It is quite intu- 
itive that in asymptotic regime the role of a priori 
knowledge should be negligible, since in principle a 
minute fraction of resources might be committed to 
prccstimate the phase and compensate for the lack 
of a priori knowledge. 

For finite N, however, a priori knowledge may 
strongly influence both the optimal precision and the 
estimation strategy itself. Since the main prospects 
for applications of quantum enhanced metrology lay 
in the regime of moderate values of N, the amount 
of a priori knowledge may play a key role in design- 
ing the optimal estimation schemes. Moreover, the 
problem is highly relevant for feedback estimation 
schemes, where in each iteration one should make 
the optimal use of the information gained in the pre- 
vious steps. This paper provides a direct way to 
find the optimal estimation schemes for an arbitrary 
form of a priori knowledge. As a result, one can also 
study the transition from global to local approaches 
by tuning the a priori probability distribution from 
uniform, p(ip) — l/2w, to peaked p{ip) ~ S(ip — ipo) 
(sec [10] which is one of very few papers approaching 
the problem of phase estimation in the intermediate 
regime) . 

Let H be an N + 1 dimensional Hilbert space, and 
let \n) G H, n = 0...N, denote an eigenbasis of 
the phase shift operator U v , U v \n) = exp(— inip)\n). 
Thinking in terms of an N photon state inside a 
Mach-Zehnder interferometer, In) denotes a state in 
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which n photons travel through the upper and N — n 
photons through the lower arm, while the upper arm 
is delayed by p with respect to the lower one. Let 
p = J2nm Pm\ n )( m \ be the probe state, which un- 
der the phase shift becomes p v — U v pU^ p . Note, 
that p may represent a decohered pure probe state, 
provided the decoherence process commutes with the 
phase shift operation, which is the case for the most 
relevant decoherence mechanisms such as photon loss 
in optical and dephasing or damping in atomic sys- 
tems [4, 5]. 

The estimation strategy is defined by a POVM 
measurement 11^, jZ^ddp = 1, 11^ > 0, where 
the measurement outcome dp is at the same time the 
estimator of the phase. Hence, Tr^II^) is the con- 
ditional probability of estimating dp provided the true 
phase is p. Taking into account the a priori probab- 
ility distribution p(tp), the average cost of the estim- 
ation strategy reads: 

C= f dpdpp(<p)Tr(p v ILf)C^, (1) 

J — 7T 

where C Vi $ is the cost for guessing dp, while the true 
value was ip. The C Vt q, — (<p — dp) 2 cost function 
is only appropriate for narrow distribution around 
p — since it docs not respect periodic conditions 
C v +2-K,(p- We choose C v .(p = 4sin 2 !£ ^ £ -, as it is the 
simplest function (with the least number of Four- 
ier components) that approximates the variance for 
narrow distributions [11], and denote its average by 

Consider an auxiliary two dimensional Hil- 
bert space Ha, where we define \ip) = (|0) + 
cxp(-^)|l))/V2. Notice that 

C^ = 4[l-Tr(\p)(p\\p)(p\)}- (2) 

Substituting the above formula to Eq. (1), making 
use of the fact that Tr(AB)Tr(CT>) = Tr[(A®C)(B® 
D)] and inserting both integrals under the trace we 
finally arrive at: 

S^p = C = 4 (1 - F) , F = Tr {RM) (3) 

where 

R= f dpp(p)\p)(p\(g> Pv , (4) 

m = [ dd5\p)(p\(g>n^. (5) 

J — 7T 

Assuming the probe state p is given, the problem 
of finding the optimal estimation strategy amounts 



to maximizing the fidelity F, over operators M which 
are of the form (5). The structure of M is ana- 
logous to the structure of a separable state, with 
an additional constraint resulting from the POVM 
completeness condition: TtaM = d^sll^ = 1. 
This observation has been employed in [12] to cast 
the problem of optimal pure state reconstruction 
into maximization of linear functional over separable 
states, which in principle may be solved numeric- 
ally using semi-definite programming. Surprisingly, 
as demonstrated below, the solution to the prob- 
lem considered in this paper is found analytically, 
without resorting to semi-definite programming, and 
just with the help of simple positive partial transpose 
(PPT) necessary criterion for separability [13]. 

Let us write the a priori distribution as p(p) = 
^Zfct-oo Pk exp(ikp) where p\ = p_ k guarantees 
p{p) G R and po = I assures normalization. Both R 
and M may be written using a block form respecting 
the structure of the tensor product Ha ® H: 

(R° M _(K M?\ (6] 

r -{rI r\)> m -\mI mi) ^ 

where R) = A (i\R\j) A , M) = A (i\M\j) A . Straight- 
forward calculation of Rj using Eq. (4) yields 

f pi\n I n / pOyi / r?l\m* 1 n 

(7) 

The completeness constraint Tr^M = 1, implies 
M ° + Ml = 1. Using the fact that R% = R\ and 
Tri?* = 1/2, we get F = \ + 23? [Tr^M?)] . 

Main result. The cost of the optimal phase estim- 
ation strategy reads 

4 (I -Kill), (8) 

where \\A\\i denotes the trace norm of a matrix. The 
optimal estimation strategy itself is given by 

N 

M = J2\ ( Pk)(p k \® \ipk)(4>k\ (9) 

fe=0 

where exp(— itpk), IV'fe) are eigenvalues and eigen- 
vectors of U = VrU r , while Vr, Ur are unitaries ap- 
pearing in the singular value decomposition (SVD) 
Rq = UbAbYr- Notice that Eq. (9) provides an 
explicit recipe for implementation of the optimal es- 
timation strategy. The optimal measurement is a 
projective measurement on basis \ipk), with pk be- 
ing the estimated phase given a measurement result 
k. 
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Figure 1. (color online) Optimal phase estimation strategies for N = 3 photon probe states |i/>) and different degrees of 
a priori knowledge Pt(tp) (black, dashed — rescaled so that f>t(0) = 1) [see Eq. (12)] for t = 20 (global regime), t — 0.2 
(intermediate regime), t = 0.02 (local regime) respectively, cpu is a phase that is estimated once a measurement 
result k is obtained, while the corresponding curve depicts p(k\ip) = | (ip^ \U V \ip) | 2 — conditional probability that 
measurement outcome k is obtained if the true phase is ip. Insets in the bottom left corners illustrate parameters at 
of the optimal probe state. In the local regime only two measurement outcomes are relevant and the optimal state is 
the NOON state — a result known from the quantum Fisher information approach. 



Proof. Let M' represent the optimal measurement 
strategy, and define M" such that off-diagonal blocks 
remain the same while the diagonal blocks are inter- 
changed: Md° = Mi 1 , Mi' 1 = Mg>. M' is of course 
positive semi-definite, but this is not guaranteed for 
M", Notice, however, that M" = WM' Ta W\ where 

is a unitary operation and Ta is the 



W = 



1 



partial transposition with respect to the A subsys- 
tem. Since M' is separable, then by the PPT cri- 
terion M' Ta > 0, and hence also M" > 0. 

F depends only on the off-diagonal blocks via 
[Tr(i2jMf )] . Therefore, M" provides the same fi- 
delity as M' . Keeping fidelity the same we construct 



M = (M' + M")/2 = 



1/2 M[° 



M'^ 1/2 



(10) 



Positive semi-definiteness of M implies that 
Mf f M{° < 1/4, hence all singular values of M[° 
are no greater than 1/2. 

Let 1$ = UrKrVI, Ml = Um^mVIj be SVDs. 
Thanks to Am < 1/2, and the fact that |Tr (J7A) | < 
TrA the following chain of inequalities holds: 



[Tr (RlM?)] < |Tr (u r ArV&U m A 



mV m 



< 



(ii) 



All inequalities are saturated for M® — U/2, where 
U = V R U R . Using the eigen decomposition U = 

Sfe=o ex P( _ * < ^fc)l'0fe)(V'fc| w e arrive at Eq. (9). ■ 

In the following examples, we choose a prior p(tp) 
that is a solution of the diffusive evolution on the 
[— 7r, 7r) interval with periodic boundary conditions 
and initial 8((p) distribution: 



\ 71=1 



cos(rup)e 



(12) 



This naturally approximates the gaussian distribu- 
tion on the interval, with t = corresponding to 
6((f) and t — oo to uniform distribution: 

For the sake of clarity, in the following examples 
we consider only the decoherence-free phase estim- 
ation, where the probe state is pure p = \ip}(tp\, 

W) = J2n=o an \ n ) ( as stressed before the method 
can be applied to more general scenarios). In this 
case, (-Ro)m = \ a n a m ex P[ — ( TO — n + l) 2 *], and the 
optimal probe state corresponds to the choice of a n 
such that the trace norm of the above matrix is max- 
imized. For t = oo (global approach), (-Ro)m = 

\ot n a n -x8 n -l,m, ll-Rolll = lEn=o' I a n 1 1 "n+l I which 
is maximized by the Berry- Wiseman (BW) states 
[11]. For general t there is no analytical formula for 
the trace norm as a function of a n and one needs to 
look for the optimal probe state numerically. 
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To facilitate the numerical search, one may imple- 
ment an efficient iterative procedure. In the i-th step 
starting with a pure probe state = \ip^)(ip^\, 
one calculates the optimal estimation strategy M^> 
- as shown above this requires only a single run 
of SVD algorithm. Having found one uses 

Eq. (3) and looks for the maximizing F — 

Tr(i?( l+1 )M«) = where 

97i w = J d<pp(<p)((p\l® Ul M {t) l®U v \<p). (13) 

Optimal is the eigenvector corresponding to 

the maximum eigenvalue of the matrix SETt^ . 

The structure of the optimal estimation strategy 
and the optimal probe state is depicted in Fig. 1 
for N = 3 and three different a priori distributions. 
Notice, how, with the increasing a priori knowledge, 
the optimal probe state evolves to the NOON state 
[1], which is the optimal solution of the QFI ap- 
proach. The periodic (2w/N) structure of the con- 
ditional probabilities p(k\ip), visible in the local re- 
gime, clearly reminds of the fact that this estimation 
strategy is useless unless the prior is highly peaked 
since otherwise there is strong ambiguity in using the 
measurement result to estimate the phase. 

It is worth mentioning, that although both the 
local approach discussed in this paper (correspond- 
ing to the limit t — > 0), and the QFI approach yield 
the same optimal probe states and the same optimal 
measurements, they in general yield different estim- 
ation precisions. While Fq is an extremely useful 
tool, it just provides a CR bound on the achievable 
estimation precision Sip > 1 / -J Fq . This bound, in 
general, cannot be saturated by an estimator based 
on the results of a single measurement. Only in the 
asymptotic limit of infinitely many repetitions of the 
experiment one may construct a max-likelihood es- 
timator which saturates the bound. The approach 
presented in this paper, on the other hand, gives an 
operationally meaningful answer for single shot es- 
timation procedure. Moreover, when t — > 0, then 
Sip — >• by the obvious fact that that the phase is 
known perfectly. Hence, for t small enough the CR 
bound is violated (this is no contradiction since our 
estimator is not locally unbiased [14]). 

Performance of the optimal estimation strategy 
for N = 10 and different probe states is depic- 
ted in Fig. 2 as a function of a priori uncertainty 

<V P rior = J f dy>4:SxnL 2 ((p/2)pt(ip). NOON states are 
optimal up to a threshold (which scales as l/N) 
above which they become useless due to 2tt/N phase 




0.05 0.10 0.15 0.20 0.25 0.30 0.35 



Figure 2. (color online) Relative reduction of the uncer- 
tainty Sip with respect to 5<p pl i OI after the optimal estim- 
ation procedure is applied to the N = 10 optimal probe 
state (black, solid); NOON state (optimal in the local re- 
gime) (red, dotted); BW state (optimal in the global re- 
gime) (green, dashed); ,, classical" state with af = (v) 
(blue, dash-dotted). 

estimation ambiguity. The optimal states clearly 
demonstrate their superiority over the BW states for 
moderate 5<^ pr i or and lose their advantage with the 
increasing prior ignorance on the phase value. 

In summary, the problem of optimal phase estim- 
ation with arbitrary a priori knowledge has been 
solved analytically, allowing to investigate the regime 
of estimation not accessible via neither the QFI ap- 
proach nor covariant measurements. Optimal meas- 
urements and estimators have been explicitly con- 
structed so that it is immediate to apply the results 
to any practical phase estimation problem. 
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